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We apply a previously developed scheme to consistently include the twist-3 distribution amplitudes 
for transversely polarized p mesons in order to evaluate, in the framework of fcr factorization, the 
helicity amplitudes for exclusive leptoproduction of a light vector meson, at leading order in Oa . We 
i-^ . compare our results with high energy experimental data for the ratios of helicity amplitudes Tii/Too 

and Toi/Too and get a good description of the data. 

QJ , PACS numbers: 13.60.Le, 12.38.Bx, 25.30.Rw 
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> ■ I. INTRODUCTION 

' Exclusive leptoproduction of vector mesons has been the subject of significant progress in the last 20 years, 

in particular, in the hard regime where a highly virtual photon exchange allows one to separate a short distance 
dominated amplitude of a hard subprocess from suitably defined hadronic objects. Experimental knowledge 
i has been gathered in a wide range of center-of-mass energies, from a few GeV at JLab to hundreds of GeV 
at the HERA collider. Following the pioneering NMC [Tl and E665 [?] experiments, the HERA collaborations 
HI and ZEUS measured p- meson electroproduction Js]. COMFASS also measured the same reaction in an 
intermediate energy range Lower energy data have been extracted at HERMES 6, "^l and JLab 

The HI and ZEUS collaborations have recently provided a complete analysis [9i, JH] of spin density matrix 
elements describing the hard exclusive productions of the p and the (j) vector mesons V in the process 



in 



X 



7*(A^)p^U(Ay)p, (1) 

which can be expressed in terms of helicity amplitudes TavA^ (A-y, Xv '■ polarizations of the virtual photon and 
the vector meson). 

The ZEUS collaboration Q has provided data for different photon virtualities —Q^, i.e. for 2 < < 160 
GeV^, 32 < W < 180 GeV ( |t| < 1 GeV^), while the HI collaboration 10] has analyzed data in the range 
2.5 < Q2 < 60 GeV^, 35 < < 180 GeV {\t\ < 3 GeV^) , where W is the center-of-mass energy of the virtual 
photon-proton system. 

The main features of the HERA data are as follows: specific scaling, and t and W dependence of the cross 
sections (features distinct from those in soft diffractive reactions) strongly support the idea that the dominant 
mechanism of the diffractive process ([T]) is the scattering of a small transverse-size, ^ 1/Q, colorless dipole on 
the proton target. This justifies the use of perturbative QCD methods for the description of the process (P). 

On the theoretical side, three main approaches have been developed. The first two, a fcT-factorization approach 
and a dipole approach, are applicable at high energy, W ^ Q ^ ^qcd- They are both related to a Regge 
inspired /cT-factorization scheme [ll| - [l7j , which basically writes the scattering amplitude in terms of two impact 
factors one, in our case, for the "f* — p transition and the other one for the nucleon to nucleon transition, with, at 
leading order, a two "Reggeized" gluon exchange in the i-channel. The Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
evolution, known at leading order (LLx) [18- 21] and ncxt-to-lcading (NLLx) order [22-25], can then be applied 
to account for a specific large energy QCD resummation. The dipole approach is based on the formulation of 
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similar ideas, not in kr but in transverse coordinate space [26|,|27|; this scheme is especially suitable to account 
for nonlinear evolution and gluon saturation effects. The third approach, valid also for W ^ Q, was initiated 
in [28^ and [1^. It is based on the coUinear QCD factorization scheme [sol l3l|: the amplitude is given as a 
convolution of quark or gluon generalized parton distributions (GPDs) in the nucleon, the p-meson distribution 
amplitude (DA), and a perturbatively calculable hard scattering amplitude. GPD evolution equations resum 
the collinear gluon effects. The DAs are subject to specific QCD evolution equations [32| - [3^ . 

Though the collinear factorization approach allows us to calculate perturbative corrections to the leading 
twist longitudinal amplitude (see [ssj for NLO), when dealing with transversely polarized vector mesons, one 
faces end-point singularity problems. Consequently, this does not allow us to study polarization effects in 
diffractive p-meson electroproduction in a model-independent way within the collinear factorization approach. 
An improved collinear approximation scheme has been proposed based on Sudakov factors 13 6l. w hich allows us 
to overcome end-point singularity problems, and has been applied to p-electroproduction [37h4C|| |. 

In this study, we consider polarization effects for reaction ^ in the high energy region, s = W"^ ^ ^ 
Aq^j^, working within the fc^-factorization approach, where one can represent the forward helicity amplitudes 
as^ 

J (k s 

where ('*'~')^''('^p)(k^, Q^) is an impact factor describing the virtual photon to p-meson transition, and 
is an unintegrated gluon density, which at the Born order is simply related to the proton-proton 
impact factor, as described below. Here, k is the transverse momentum of the channel exchanged gluons. 

The impact factor (-^7)-^p(-*'p)(k^, Q^) vanishes at k — >■ 0, which guarantees the convergence of the integral 
in Eq.([21) on the lower limit^. In fact, ^-^^^p^ ^ k^/Q^ at |k| ^ Q, which allows us to express the longitudinal 
amplitude in the collinear limit in terms of the usual gluonic parton distribution function. 



2. f dk, 
xg[x,Q ) = / -^T[x,k) 



In the case of the transverse amplitude, the situation in the collinear limit is different due to the above- 
mentioned end-point singularities. At |k| <^ Q the transverse impact factor $^«_>.py ~ (k^/Q^) In (Q^/k^) and 
the transverse amplitude cannot be expressed in terms of the gluonic parton distribution function. Nevertheless, 
both longitudinal and transverse amplitudes can be calculated within the fc^-factorization description For 
the transverse amplitude the end-point singularities are naturally regularized by the transverse momenta of the 
t— channel gluons [43 45]. 

Moreover, at large photon virtuality providing the hard scale, the 7* — p impact factors can be calculated in 
a model-independent way using QCD twist expansion in the region k^ 3> Aq,^^,. Such a calculation involves 

the p-meson DAs as nonperturbative inputs. The principle point here is that the region k^ ^ ^qcd gives the 
dominant contribution to the amplitudes in the integral ([2]) . Below we introduce an explicit cutoff for transverse 
momenta of i— channel gluons to clarify this point. The calculation of the impact factors for ^^^^p^^, <I>^j,_j.p^ 
is standard at the twist-2 level (461], while $^^_j.py was only recently computed [44^7 45] (for the forward case 
t = tmin), up to twist-3, including two- and three-body correlators, which contribute here on an equal footing. 

In our study we use results [4^, for the ^-y^-s-pi, Etnd $^^_j.py impact factors and a phenomenological 
model (47! for the proton-proton impact factor. This model involves a single energy scale parameter M, and 
is equivalent in our calculation to a specific assumption for the k dependence shape of the unintegrated gluon 
distribution J-"(x, k). 

Our approach is close in spirit to the calculations performed in Refs. [48l-l52| within the dipole approach, 
where amplitudes are related to the light-cone wave functions (j){z, r). Collinear DAs are integrals of light-cone 
wave functions in momentum space over the relative transverse momentum conjugated to r. The light-cone wave 
functions are complicated objects, and in practice, their dependence on the longitudinal momentum fraction 



^ We use boldface letters for Euclidean two-dimensional transverse vectors. 

^ This property of the impact factors is universal and related to the gauge invariance [4ll. |42|| . It is a consequence of QCD gauge 
invariance and is in accordance with the general statement of the Kinoshita, Lee and Nauenberg theorem which guarantees the 
infrared finiteness of amplitudes in the case of the scattering of colorless objects. 
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z and transverse coordinate separation r variables, which describe the gg dipole, should be modeled. Our 
main point here is that one can assume the dominant physical mechanism for production of both longitudinal 
and transversely polarized mesons to be the scattering of small transverse-size quark-antiquark and quark- 
antiquark-gluon colorless states on the target. This allows one to calculate corresponding helicity amplitudes 
in a model-independent way, using the natural light-cone QCD language - twist-2 and twist-3 DAs. 

This paper is organized as follows. In Sec. |lTl we recall some results of Refs. [45| and [HI] about the impact 
factors and distribution amplitudes for the meson, which we need for computing helicity amplitudes. In 
Sec. mil we describe first the model for the proton-proton impact factor. Then we compare with HERA data 
the ratios of helicity amplitudes Tn/Too, calculated both in the Wandzura-Wilczek (WW) approximation and 
with the account of the genuine contribution. We also compare our predictions for the Tqi/Toq ratio with the 
data of HERA. We obtain a good description of these two ratios. For both observables, we discuss the effect 
of the energy scale for the proton-proton impact factor, as well as the sensitivity to the infrared region of 
t- channel gluon momenta. 



II. IMPACT FACTORS 7* (At) ^ p(Ap) 



A. Impact factor representation 



In the impact factor representation at the Born order, the amplitude of the exclusive process ^*(\^)N — )■ 
p(Ap) N reads 

r,^,^(r;Q,A/) = zs / ^^^^$^^-^(k, r; A/2)$7*(a. W(a,) ^. q2) ^ (3) 
J (^T^) k (k — r)^ 

as illustrated in Fig. [TJ Note that in this representation, no skewness effect is taken into account. The 
7*(A-y) — )■ p(Ap) impact factor ^'^ (^y)'^p'^^p^ is defined through the discontinuity of the S matrix element 
for 7*(A^;g)g(fc) g{r - k)p{Xp;pp) as 

^,-(A,)^p(A,) ^l_J ^j^.^^^ , (4) 

where k = (k + q)^ . In Eqs.Q and Q the momenta q and pp are parametrized via Sudakov decompositions in 
terms of two lightlike vectors pi and p2 such that 2pi.p2 — s, as 

q = Pi P2 and p^pi^ p2 + rj_ , (5) 

s s 

where = —q^ » Aq^^, is the virtuality of the photon, which justifies the use of perturbation theory, and 
rUp is the mass of the p meson. The impact in Eq.Q cannot be computed within perturbation theory, 

and we will use a model described in Sec lIII Al Note that due to QCD gauge invariance, both impact factors 
should vanish when either one of the transverse momenta of the channel exchanged gluon goes to zero, k — > 
or r-k — )■ . 




Figure 1: Impact factor representation of the 7* A — >• pA' scattering amplitude. 
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The computation of the 7* — )■ p impact factor is performed within cohinear factorization of QCD. The 
dominant contribution corresponds to the 7^ — > transition (twist 2), while the other transitions are power 
suppressed. The 7^ — pL and 7^ — )■ pL impact factors were computed a long time ago [46j. while a consistent 
treatment of the twist-3 7^ — > pT impact factor has been performed only recently in Ref. |45| . It is based on 
the light-cone coUinear factorization (LCCF) beyond the leading twist, applied to the amplitudes 7*(A^)(7(fc) 
g{r—k)p{Xp), symbolically illustrated in Fig.[2l Each of these scattering amplitudes is the sum of the convolution 




Figure 2: Two- and three-parton correlators attached to a hard scattering amplitude in the specific case of the 7* ^ p 
impact factor, where vertical lines are hard t— channel gluons in the color singlet state. 



of a hard part (denoted by H and Hp for two- and three-body contributions, respectively) that corresponds to 
the transition of the virtual photon into the constituents of the p meson and their interactions with ofF-shcll 
gluons of the t channel, and a soft part (denoted by $ and As the photon is highly virtual, this convolution 
reduces to a factorized form, expressed as a convolution in the longitudinal momentum of hard parts in collinear 
kinematics with DAs. 



B. Distribution amplitudes in the LCCF parametrization 



Our accuracy is limited to dominant contributions both for 7^ — p^ (twist-2) and 7^, — px (twist-3) 
transitions; therefore, only leading terms of the expansion in l/Q in both amplitudes are kept. Hence, only two- 
body (quark-antiquark) and three-body (quark-antiquark gluon) nonlocal operators are involved. Correlators 
and distribution amplitudes depend on a factorization scale p. We have to take into account this dependence 
to compare our results with experimental data. 

The seven chiral-even p-meson DAs up to twist 3 are defined by matrix elements of nonlocal light-cone 
operators. Let us introduce p and n, two light-cone vectors such that pp ^ p at twist 3 and p ■ n = 1. The 
polarization of the out-going p meson is denoted bye*. 

The two-body correlators are parametrized as'* [45| 



(6) 
(7) 
(8) 
(9) 



{piPp)\^iz)jpm\o) - rnpfp dye^yp-^[My;t^^W-n)p^. + My;f^^VTp] 

Jo 

{p{Pp)\i^{zh5lf.m\0) = impfpep^pse*Ty^n' ( dy e^yP-'My; p^) 



p{Pp) 



p{Pp 



= mpfpPpe*T^ / dye'yP-'^l{y;p^) 



i-mpfppp£a\!3se*r^p^ ■n!' I dy e^y'^-'' ipl^{y; p^) , 



where y and y — 1 ~ y are, respectively, the momentum fractions of the quark and the antiquark, while the 



^ The chiral-odd twist-2 DA for the transversely polarized p meson does not contribute to the process considered at the accuracy 

discussed here. This is also true in the approach based on collinear factorization of generalized parton distributions [s^ . |55| . 
* In the approximation where the mass of the quarks is neglected with respect to the mass of the p meson. 
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three-body correlators are expanded as 

/■I pV2 

{p{Pp)\i'{zih^gAl{z2)m\0)^rn,fl{f,')p^e*rJ dy^ dy.e^y^P-^^+^^y^-y^^P-^^ B{y,,y,; fi^) (10) 

Jo Jo 

fi fV2 

Jo Jo 

where yi, ?/2, and yg ~ y2 ^ J/i are, respectively, the momentuni fractions of the quark, the antiquark, and 

the gluon. We used the standard notation dp = ■^{dp — dp) . Normalizations and symmetry properties of 
the distribution amplitudes are recalled in Appendix \^ and they are used to simplify helicity amplitudes in 
the forthcoming computations. For later use, we define the following combinations of three-body distribution 
amplitudes: 



M{y,,y2;p^) - (Hf^) Biy^^yr, p') - Ciip') Diyi,y2; p') , 
5(2/1, ya;^') = C3p{^^^)Biyl,y2;^i')+C3p{^'^)Diyl,y2;^i') 

where CspiP^) a-nd Cs^Cm^) dimensionless coupling constants: 

.V,2,_fM^^ ,A,2,_fiM) 



fp ^ fp 

C. Reduction of DAs to a minimal set ifi, B, D 



(12) 



(13) 



DAs are not independent; they are linked by linear differential relations derived from equations of motion 
and n independency [H,!!!]. The solutions for (pp{y) = {(p3,(pA,'fi ,f'A} are the sum of the solutions in the 
so-called WW approximation and of genuine solutions: 



^P(y)=^^^(y) + ^r(2/)- 



(14) 



The WW approximation consists in neglecting the contribution from three-body operators by taking 
B{yi,y2]iJ?) = D{yi,y2;p^) = 0. tp^"^ {y) are then only functions of tfi{y): 



WW / 2\ 
WW I 2\ 

^27(2/; M^) 
^1^^(2/;m^) 



u 
u 



du 



2M 



V 



u 



u 



Jo U jy U 

yf\u^l^+yf\u^^^-'^'^ 

Jo U jy 



U J y U 



gen/ 2\ 



qen f ^\ 



du 



du 



A{u;n^) 



du 



A{u; p^ 



where A(u] jj?) has the compact form 



A(u;n'^ 







/ dy2 


— du 


Jo 


J/2 - M 



f 1 

M{y2,u;p^)+ / dy2 M{u,y2;p^ 

Ju 2/2 - U 



(15) 
(16) 
(17) 
(18) 



Genuine solutions only depend on {B{yi,y2; p ), D{yi,y2; p )} or the combinations 
{'5(2/1, 2/2; M^),M(yi, 2/2; M^)}: 



(19) 
(20) 

(21) 
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and it obeys the conditions 

[ du A{u; fi'^) ^ and / du u A{u; fi'^) ^ , (22) 

JO ^0 

coming, respectively, from the constraints 



^r(2/;/.^)dy-0 and / {y - y) ip^ iv; ^^') dy = . (23) 
Jo 

Equations (|19p and I^U^ determine the expressions of if\^{y\ fi^) and (p^j^^{y; /i^) as 

gen, 2n r, gen, 2n ^ ^(yi , ^2 ; M^) + -^^^(^1 , 2/2 ; A^^) 

nriy^f^) = / ^^"'/'a Wn)~- dyi dy2 , (24) 







2 Jo Jy y2 - yi 



gen, 2^ T , gen, 2n 1 /"^ , '5'(?/l , ^2 ; M^) " M {yi, y2; 







2 Jo Jy V2 - yi 



The correspondence between our set of DAs and the one defined in Ref. [5j| is achieved through the following 
dictionary derived in Ref. [i^. It reads, for the two-body vector DAs, 



and for the axial DA, 



vi{y)^h'^y), V3iy) = 9±Hy), (26) 



(2r) 

For the three-body DAs, the identification is 

uf \ ^(2/1,1-2/2) , ^ A{yi,l-y2) , , 

B{yi, 2/2) = and £1(^1, j/2) = • (28) 

2/2 - 2/1 2/2 - 2/1 

Explicit forms for ipi, B, and D are obtained with the help of the results of Ref. [53*1 obtained within the QCD 
sum rules approach. The first terms of the expansion in the momentum fractions of the three independent DAs 
thus have the form 

My,^^^) = 6ro(l-f a2(/i')^(5(y-y)2-l)), (29) 

5(2/1, 2/2; M^) = -5040yij/2(j/i - 2/2)(2/2 - 2/1) , (30) 

i^(yi,2/2;M') = -360yiy2(2/2-2/i)(l + ^i^ (7(2/2 - 2/i) - 3)) . (31) 

The dependences on the renormalization scale of the coupling constants 02, uj^^ qj, Ca^, and are given 
in Ref. [s^. In Appendix B we present both the evolution equations and the values of these constants at 
fx^ = \ GeV^ used in our analysis, as well as the dependence on jj? of the DAs. 

D. Impact factors 

In the Sudakov basis, the longitudinal and transverse polarizations of the photon are^ 

<L = 7^K + — P2), e± = 4s(0'Tl,-*,0) . (32) 



In Ref. [45| we took = (0, 1, iti, 0), which we change here for consistency with the usual experimental conventions [56 
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For t = tjnin the same parametrization will be used for the p-meson polarization with — >■ —m^ and Q — >■ irip. 
We introduce the notations: a — ^ and B = 2'Kas-^fp. The impact factor 7^ — ^ has been computed up 
to twist 2; the next term of the expansion is of twist 4. It reads [i^ 



^,,^p, (k, r, Q; ^ 2QB^^ / ^ ^ ^' 



with 



i^p(y,k,r,Q) 



1 



1 



1 



(33) 



(34) 



Here a and 6 are indices of color and Nc is the number of colors. 

The first non- vanishing term of the power expansion of the ^ pl impact factor has been calculated in 
Ref. [4^. It corresponds to the twist 3 and, in the limit t — >• imin, reads 



*7T^Pi,(k,r,Q;p^) = J dy (y - y)<^i(y; p^) e • Qp(?/, k, r, Q) 



with 



Qp(y,k,r,g) = 



yr 



yr 



k — yr 



k-yr 



(yry + yyQ^ {yr^ + yyQ'^ (k - yr)2 + yyQ2 - yr)2 + yyQ2 



and where 



(35) 



(36) 



n/2 



(Tl,-*) 



The impact factor for — )■ pT with the exchanged momentum r = is [45| 



2^ _ (ei-ep ^BrUpSab 



a{a + 2yy) ^ 2n , 2m 

- dy— — ■ — —:^[{y-y)ip-i^(y;p) + ipj^(y;p) 

yy(a + yy) 



yiyia 



1 /■a2 
c^2/2 / dyi' 

Jo " + yiyi 

dy2 / rfyi 
Jo 



2 - A^JC; 



1 



2 + jV,/CF 

A^c (2/2 -yi)y2 



a(yi + 2/2) + yiy2 C f yia + yi{y2 - yi) 
yi / (2 - Ar,/Ci.)yia 



M{yi,y2;p^) 



a + yiyi \a{yi + ^2) + yiy2 



- 2 



1 



yi 



ayi + (2/2 - yi)y2 



S{yi,y2;p'^ 



(37) 



(38) 



where Cp = ^^f^p- is the quadratic Casimir of the fundamental representation. One can readily check from 
Eqs. (|33|) and ([35]) that $-y^_j.pj- (k, r, Q; /x2) and ^^^-:,p^ik,T,Q\ p^) vanish when k — > or k — > r. Similarly, 
one can see from Eq. ( l38| that $-y^_).pj, (a, Q; p?) vanishes in the limit a — >■ 0. 



III. HELICITY AMPLITUDES 

A. A phenomenological model for the proton-proton impact factor 

To compute ratios of helicity amplitudes, we need a model for the proton-proton impact factor. A simple 
phenomenological model was provided for hadron-hadron scattering in Ref. [47| , of the form 



1 



1 



M2 + (|)2 M2 + (k-§)2 



(39) 



A and M are free parameters that corresponds to the soft scale of the proton- proton impact factor. We discuss 
later the value of M . The parameter A has no practical importance for our study since the observables we will 
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be interested in involve ratios of two scattering amplitudes, which are insensitive to A. Note that this impact 
factor indeed vanishes when k — or r — k — > in a minimal way. 

This simple model can be interpreted by assuming that, inside the proton, there exist some typical color- 
dipole configurations (onia) which will couple to the p— meson impact factor through a two-gluon exchange. 
This impact factor has the same form as a 7* — >■ 7* impact factor, with a scale which governs the typical 
transverse momentum. Such a model was the basis of the dipole approach of high energy scattering [57| and 
used successfully for describing deep inelastic scattering at small x [58| . 

Returning to the discussion given in the Introduction, one can reformulate this model for the proton impact 
factor at r = into a simple assumption about the form of k dependence of unintegrated gluon distribution, 

k^ 

T{x,k) 5 . (40) 



B. Helicity amplitudes Tn and Too at t = tmin - Comparison with HERA data 



Let us start from the impact factor representation at the Born order given by Eq. ^ . The helicity amplitude 
Too is, from Eqs. (0 iMl [Ml [SI) , 



To, 



isCp^AB 
(27r)Q5 



dyviiy,^'^) / da 



1 / 1 



2 \ i?^ a + R"^ J a -\- yy 



(41) 



In order to obtain the WW contribution to the Tn amplitude, three integrals should be performed. One, 
over a, is related to the transverse momentum of t— channel gluons, and two, over y and u, are related to DAs; 
see Eqs. (|T7| and (fT8)) . It is useful to interchange the order of integrals over a, y, and u in order to fix a specific 
model for DAs at the last step when performing the u integration. This leads to 



WW 



where R 



rpVV 
^11 



2 _ 



isCF{e^.e;)mp2AB r\ (^i(u;^2) 



(2^)g6 



du- 



1/1 1 \ a{a + 2yy) 



(42) 



We also introduced a cutoff i?f = ^ on the integral over a , where A is the cutoff on |k|. This 



cutoff allows us to see how soft gluons contribute to the amplitude. 
The genuine contribution is 



.gen isCF{e^.e*p)mpAB 
11 — 



(27r)Q6 



da 



1 

i?2 



1 



i?2 



dy2 / dyi 



y2 



dy2 I dyi 



a + yiyi 
2 + NJCf 

m 



2 - NjCf 



Nr. 



a{y\ + yi) + yiy2 Cp y2a + yiiy2 ~ yi) 



a(^a ^ 2yy) _ gen/ 2\ , gen > 2\i 

dy—, — ■ — -2 (y - ywiT M + vat y; ) 
yy{a + yyr 

1 

M(yi,2/2;/x") 



2/1 



Ne {y2 - y\)y2 



Cp 



(2 - Ne/CF)yia 
a + yiyi \ a{yi + ^2) + 2/12/2 



5'(yi,2/2;/^^) 



2/1 ayi + (2/2 - 2/1)2/2 



(43) 



For convenience, we define Ii{y\R^,R\), ^2(2/1: 2/2; ^1), and I3{yi,y2; R^ , Rl) as the integrands after inte- 
gration over a: 



h{y;R\Rl) 
l2{yi,y2;R\Rl) 
h{yuy2;R\Rl) 



da 



1 



a{a + 2yy) 



da 



1 

'r^~ 
1 

i?2 a + i?2 J a + yiyi 



J yy{a + yy)^ 



2/12/1" 



2 - Ne/Cf 



Nr. 



a(2/i + 2/2) + 2/12/2 2/2" + 2/1(2/2 - 2/1) 



da ■{ — - — 



1 



i?2 a + i?2 



2 + Ne/Cp 



2/1 



(2 - Ne/CF)yia 
a + 2/12/1 V"(2/i + 2/2) + 2/12/2 



2/1 



(2/2 - 2/1)2/2 



1 



Cf 2/1 Oiyi + (2/2 - 2/1)2/2 



(44) 
,(45) 

(46) 
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such that can be expressed as (removing the variables and for simphcity) 

{27t)Q' 



rr^gen IS C pie-y ■e*p)mpAB ( 

Til = T^T^^Tvi i~ / dy Ii(y)[[y - y)iplj, (y; n ) + ip\j,{y; )\ 







(•y2 /■! /■a2 

+ / dy2 I dyi/2(yi,2/2)M(yi,y2;Ai^) - / ^2/2/ ^2/1/3(2/1, y2)5'(yi, 2/2; M^) f- , (47) 

^0 Jo 



which, using the symmetry property 5(2/1, 2/2; A*^) = —M{y2,yi;fJ.^), turns into 



1 



T{i = ^^^^Qg \- I dyli{y)[{y-y)iplj,{y]^i^) + ip\j,{y]^i^)] 







1 /■1/2 

dy2 f dyi{l2iyi.,y2) + l3iy2,yi))M{yi,y2;n^)} (48) 
Jo 



with 



^2(2/1,2/2) +^3(2/2,2/1)= 2-— da- 



CfJJrI -R^(a + i?2)(a(2/i + 2/2) + 2/12/2) V" + 2/22/2 a + yiyi 

H / c?Qf — 3- ^-7 —7 H / da — — — — ■ . (49) 

Cf JrI R^{a + R^){a + yiyi){ay2+yi[y2-yi)) y2 JrI i?^ (a + i?2)(a + 2/22/2) 

Combining the results gj) with g2]) and (08]), the ratios T^^ /Tqq and Tf^/Too read 

1 , 2, f\ r , a + 2xv{l-xv) f 1 1 
dv ipiyv] fi j dx da- ' 



M Al - "^P ^o ^ Jflf a(« + x^(l-x^))Hfi^ " + ^^ 

dy<pi[y,n ) 







a (a + 2/2/) a + 



where we took into account that £'y£*p — —£7 • — —1, and 



1 pi <.y2 

9<="('„-,,2^ 1 ,„9eri/ 2ni_ / J„„ / rl.,,(T^(o,. n,„\j.T„(^,„ n-,A\T\/r(n,. ,,„-,,2^ 



^ , dy Ii{y)[{y - y)^pl^ {y; n ) + ip^^^ {y; H )] - dy2 ^2/1 (-^2(2/1, 2/2) + ^3(2/2, 2/i))^(2/i, 2/2; ) 
^ 11 _[^Jo Jo Jo 

Too 2g f\ ^ 2. r ^ ^ f ^ 1 \ a 

/ dyipi(y;fj,) da— I — — ^ — ■ — = 

Jo JRj a \R a + R^J a + yy 

(51) 

The integration is performed analytically over a and numerically over the variables left as, for example, y for 
Too, X, V for T]^^ , and 2/1, 2/2 for Tf"^ . The measured ratio Tn/Too is conventionally defined [Hg'I to have the 
opposite sign with respect to Eqs. (1501) and (1511) . in order to ensure the usual matrix summation in the definition 
of the density matrix. 

In Fig. [3] we show different contributions to the ratio Th/Tqo as a function of and for typical values of 
nonperturbative parameters M and A. Unless specified, we take as a factorization scale fi — Q. Note that the 
factorization scale only appears in the ratio of the amplitudes through the DAs and the coupling constants. We 
see that the WW contribution dominates over the genuine one. For illustration, we also show this ratio using 
the asymptotic ipl'^ — 6yy DA, which corresponds to /i — > 00, which is identical to the WW contribution since 
the genuine twist-3 contribution vanishes in this limit. The small difference between this last curve and the 
solid one shows a weak dependence of this ratio on the factorization scale ^. 

The two parameters A and M have different physical meanings. M is the typical nonperturbative hadronic 
scale, while A is the minimal virtuality of gluons, which should be bigger than Aqcd for consistency of our 
perturbative approach. From Fig. |4] (left panel), we see that our predictions are stable for M in the range 1-2 
GeV. The data, when compared with our model, with n = Q, favor a value of M of the order of 1-2 GeV but 
exclude a very small value around Aqcd- From Fig. [J] (right panel), we see that for A around Aqcd, our results 
are very close to the experimental data and rather stable, whereas for A = 1 GcV, i.e. significantly larger than 
Aqcd — 220 MeV in the MS scheme, they notably deviate from the data. Let us stress that the fact that our 
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+ 

0.2 - 



^ -r7---7-^-rT^- , -rTn- , ~r.- , - , 

5 10 15 20 25 

Figure 3: The WW contribution t}^^ /Too in green (dash-dotted line), the genuine contribution Tf^"/Too in red (dashed 
line), and the sum of the two contributions in blue (solid line), at M = 1 GeV and A = GeV, as functions of the 
virtuality of the photon. The brown (long-dashed) curve is the contribution based on the asymptotic DA of the p 
meson, ipi{y,fi^ = oo) = ipi^i^y) = 6j/(l — y). Our results are compared with the experimental data from HI The 
experimental errors are taken to be the quadratic sum of statistical and systematical errors. 




; t ~ -4-.-. 

0.2 0.2 



5 10 15 20 25 5 10 15 20 25 

Figure 4: Predictions for the ratio Tii/Too as a function of Q^, compared to the experimental data from HI [l^. The 
experimental errors are taken to be the quadratic sum of statistical and systematical errors. Left panel: Fixed A = 
GeV cutoff and various values for M. Right panel: Fixed scale Af = 1 GeV, and various values of the cutoff A. 



estimate provides the correct sign for the ratio Tqi/Too when compared to HI data is a nontrivial success of our 
approach. 

In Fig. Owe show the result of our calculations for the spin density matrix element Tqq as a function of and 
for typical values of the nonperturbative parameter M and for A = GeV. This observable allows a comparison 
of our prediction with the whole set of HERA data^. Note that our amplitudes are evaluated at t = tmin while 
experimental data are integrated over some t range but dominated by very small values of t. Ai t — t^in, Vqq 
only depends on the ratio xu — |Tii|/|Too| through 

'^00 = ^ , 2 ' (52) 



We predict ratios of amplitudes, while ZEUS made available the spin density matrix elements; HI extracted both spin density 
matrix elements and ratios of amplitudes. 
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where e is the photon polarization parameter e ~ (1 — y)/(l — y + 2/^/2). For HI (e) = 0.98 and for ZEUS 
(e) = 0.996. 

mm 7 ' 00 

Experimental data are dominated by \t — tmin\ < 0.4 GeV^, for which the only significant amplitudes are 
I Too I > \Tii\ > I Toil. The exact relation reads 



For t 7^ tmim i"oo Only slightly depends on the s-channel helicity violating amplitudes Tqi, Tio, and Ti-i. 



= ^ + ^01 (53) 

•^11 £ ~\~ ^01 "^1 1 "^lO 

where Xij — |Tij|/|Too|. The expression (j53p simplifies when one neglects the contributions of xi-i and xio, 

x'^ ■ 

and expands the remaining expression up to the first power of . In this way, we get 



„04 ^ 1 / ^ I -^Ol-^ll 



(54) 



Based on HI and ZEUS data, the effect of the second term in the parenthese of Eq. (|54p is below 1 %, so up to 
this accuracy the expression ([54|) is equivalent to the formula (|52|) . 



A = GeV 

M = 0.5 GeV 
M = 1 GeV 
M = 2GeV 



10 20 30 40 50 



Figure 5: The spin density matrix element rJJo as a function of for M — 0.5 GeV in green (dash-dotted line), M — 1 
GeV in blue (solid line), and M — 2 GeV in red (dashed line), and for A = GeV. Our results are compared with the 
experimental data from ZEUS ^ and HI [l^]. The experimental errors are taken to be the quadratic sum of statistical 
and systematical errors. 



C. Helicity amplitudes Too and Toi for t ^ tmin 

The HI data show that the spin-flip amplitude Toi is nonzero, showing an explicit s— channel helicity violation. 
Besides, this amplitude vanishes when the squared momentum exchanged by the proton t = — is zero. We 
start with the generalization of Eq. (|41|) for t ^ tmin ■ 



Similarly, 



isCp^QAB r J... 2n /■ (k- r/2)2 - (r/2)2 



(27r)2(M2 + (r/2)2)7o ^ ' 7 k2(k - r)2 (k-r/2)2 + M2 

111 1 



{yr)^ + yyQ^ {yr)^+yyQ^ (k-yr)2 + yyQ2 (k - yr)2 -|- ( 



)2 



zsCf2QAB f\ ^ f d^k (k-r/2)2-(r/2)2 

Tqi = /o \9rn.r9 , r / dy [y - y) ipi{y; n ) 



(27r)2(Af2 + (r/2)2) Jo " ^'^""^ ^ k2(k - r)2 (k-r/2)2-HM2 

yr-e yr^e (k— yr)-e (k— yr)-e 



(yr)2-KyyQ2 (yr)^ + yy (k - y r)2 yy Q2 (k _ yj.)2 ^ g2 



(55) 



(56) 
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The kT integrations in Eqs. (|55l I56p are performed in two different ways, presented in detail in Appendix [Cl 

• the integrations over fcy are performed, without any infrared cutoff, partially analytically through a residue 
method (see Appendix I C 1|) . 

• the integrations over fc^ are performed, with an infrared cutoff, fully numerically through triangulation 
coordinates centered at the pole of the two channel gluons (see Appendix IC 2\ . 

Figure ini shows the dependence of the ratio Tqi/Too on the choice of the factorization scale ji for M — 1 GeV 
and A = GeV. For completeness, we also show the predictions based on the asymptotic DAs. We see that 
for factorization scales around = our results are rather insensitive to its values. Nevertheless, the ratio 
Tqi/Tqo seems to be more sensitive to this scale than the ratio Th/Toq. 




0.0 0.5 1.0 1.5 2.0 ' 0.0 0.5 1.0 1.5 2.0 



Figure 6: Predictions for the ratios Tbi/Tbo as a function of \t\ for M = 1 GeV and A = GeV, for different values of the 
factorization scale /i^, compared with HI data the blue (solid) line is for (^i(j/,^^ = Q^), the green (dotted) line is 
for ipi{y,ii^ = 2Q^), the brown (dashed) line is for (^i(j/,/i^ = and the red (dashed) line is for ipi{y,n^ — oo) = 

ft^iv) = 6y(l — y). The experimental errors are taken to be the quadratic sum of statistical and systematical errors. 
Left panel: = 3.3 GeV^ Right panel: = 8.6 GeV^ 

Our predictions are based on perturbative QCD and therefore, at small t, can only lead to a powerlike or 
logarithmic t dependence. On the other hand, it is well known that data exhibit an exponentially falling t 
distribution, as has been seen both by the HI [lO^ and ZEUS Q collaborations. Therefore, we have studied the 
effect of multiplying our predictions for the amplitudes by a factor e~^^ l*-tminl/2^ where bi {i = L, T) corresponds 
to p electroproduction from 7^ or 7^, . HI measured values of 6l and 6l — 6t [lOj. The measured values for 
the latter are - br = -0.03 ± 0.271°;^? GeV-^ (for (Q^) ^ 3.3 GeV^) and bL ~ br = -0.65 ± O.Wt^ ji 
GeV~^ (for (Q^) — 8.6 GeV^). Here we present our results in Fig. [71 One can see in the right panel of Fig. [7] 
that the precision of the data for the Tgi/Too ratio does not permit us to discriminate between a zero value 
for the difference of the transverse and the longitudinal slope parameters, 6l — br, and a nonzero value of this 
difference, as measured by HI at higher values of . 

Thus our estimate provides the correct sign and order of magnitude for the ratio Tqi/Too when compared to 
HI data for M of the order of 1 GeV in the whole range of {-t) < 1.08 GeV^. 

For completeness, as we did for the ratio Th/Tqq, we also display in Fig. |S]the effect of varying the cutoff A 
on kr for the ratio Toi/Too- Again, the prediction does not change significantly when A is around Aqcd- One 
obtains the same kind of values for Al and A when comparing with the data for the two ratios Th/Tqq and 
Toi/Too- However, due to a lack of precision of the data for the ratio Toi/Too, the parameters M and A are 
mainly constrained by the ratio Th/Toq. 

IV. CONCLUSION 

We have evaluated the ratios Tn/Too and Toi/Too within kr factorization, which we compared with recent 
HI data. We got fairly good agreement with reasonable values of the two nonperturbative parameters A 
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Figure 7: Predictions for the ratio Toi/Too as a function of \t\ for A = GeV, for various values of M, compared with HI 
data The experimental errors are taken to be the quadratic sum of statistical and systematical errors. Left panel: 
= 3.3 GeV^ Right panel: = 8.6 GeV^ 



Too 



0.4 - 



0.3 



0.2 



0.1 



M = 1 GeV, = 3.3 GeV^ 

A = GeV 

A = 0.2 GeV 

A = 0.4 GeV 



0.0 



0.5 



1.0 



1.5 



2.0 




Figure 8: Predictions for the ratios Toi/Too as a function of |t| for M = 1 GeV, for different values of A, compared with 
HI data [13]: the blue (solid) line is for A = GeV, the red (dashed) line is for A — 0.2 GeV, and the green (dash-dotted) 
line is for A = 0.4 GeV. The experimental errors are taken to be the quadratic sum of statistical and systematical errors. 
Left panel: Q"^ = 3.3 GeV^ Right panel: Q"^ = 8.6 GeV^ 



and M involved in our description. In particular, we found rather weak sensitivity of the obtained helicity 
amplitudes to the region of small i-channel gluon momenta (small values of A). This justifies our approach 
to the (-'^T-)~*'^('^p)(k^, Q^) impact factor, based on the assumption about the dominance of the scattering of 
small transverse-size quark-antiquark and quark-antiquark-gluon colorless states on the target, which includes 
its calculation using QCD twist expansion in terms of meson DAs. 

Our calculation of the transverse amplitude involves both WW and genuine twist-3 contributions. It turns 
out that with the input for coupling constants in DAs of p mesons, determined from the QCD sum rules [s^ , 
the WW contribution strongly dominates these two observables. Besides, we found rather weak dependence of 
the WW contribution on the shape of meson twist-2 DAs. In other words, we found that the Tn/Too amplitude 
ratio turned out to be not very sensitive to the physics encoded in the meson DAs. This opens a possibility 
to constrain the other important quantity which enters the /cT-factorization formalism, the unintegrated gluon 
distribution J-{x, k). In the present study we use the very simple ansatz (|40|) for the k shape of F{x, k). Having 
rather precise data for the Tn/Too amplitude ratio, it would be very interesting to test different approaches to 
unintegrated gluon distribution present in the literature, calculating this ratio in fc^ factorization. 

The lack of precision for the present data for the ratio Toi/^oo does not allow us to fix precisely the non- 
perturbative parameters in our description. Nevertheless, their values [M being of the order of 1-2 GeV and A 
being of the order of 0-0.2 GeV) extracted from the Tn/Too ratio analysis and used for the prediction of the 
2oi/2oo ratio do not contradict data. 
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Other scattering amplitude ratios have also been measured and should be confronted by a /cy-factorization 
approach. This requires nontrivial analytical calculations for t ^ tmm of the twist-3 amplitudes (which was 
not needed for the ratio Tqi/Tqo) which is a hard task, since it involves, in particular, the computation of the 
72 — > pt impact factor. This deserves a separate study. 

Data also exist for (p leptoproduction. In this case quark- mass effects should be taken into account, in 
particular, because this allows the transversely polarized 4> to couple through its chiral-odd twist-2 DA. The 
fact that the ratio Th/Tqq is not the same (after trivial mass rescaling) for p and (f> mesons points to the 
importance of this effect. This is also beyond the scope of our present study, but may open an interesting way 
for accessing chiral-odd DAs. 

The BFKL resummation effects have not been included. Although they were expected [1^, [s^ to be rather 
dramatic at the level of scattering amplitudes, as was confirmed by a leading logarithmic x analysis of HERA 
data [60l - [6^ . we expect that for the ratios of amplitudes considered here, they should be rather moderate. Also, 
the next-to-leading order effects - both on the evolution and on the impact factor - should be studied. 

On the experimental side, the future Electron-Ion Collider with a high center-of-mass energy and high lu- 
minosities, as well as the International Linear Collider, will hopefully open the opportunity to study in more 
detail the hard diffractive production of mesons [63l - [70j . 
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Appendix A: Distribution amplitudes 



dy 'fiiy; iJ.'^) ^ I , I dy ¥'3(2/; M^) = 1 , / dy {y - y)ip Aiy; t^'^) ^ \ ■ (Al) 



Distribution amplitudes are normalized as 

»i /.I /.I 

d.ii (11 — Ti^in A (11: 11^^ — ^ 

The following symmetry properties are derived from C-parity analysis: 

B{yi,y2;f^'^) = -B{y2,yi;p'^) , D{yi,y2; p.'^) ^ D{y2,yi; p"^) . (A2) 
Then, the functions M and S defined in Eq. satisfy 

M{yi,y2;p^) = -S{y2,yi;fJ-^)- (A3) 



Appendix B: Evolutions of DAs as a function of Q 

The parameters entering the DAs at p = I GeV and their evolution equations are given in Ref. ^53,]. In 
Table U we recall their values for the p meson at /ip = 1 GeV. 
For 02, the evolution equation is 

a2{p') = a2{pl)L{py^-/'"> (Bl) 

with 



asipl) l + !fa,ipl)Hpypl) 
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0.52 




-2.1 


V 
'^[O.ll 


28/3 


0.2, p 


0.18 ± 0.10 


rA 


0.5 - 0.6 10-2 GeV^ 


■m-php 


0.2 10^2 GeV^ 


'-,Sp 


0.032 


(sp 


0.013 



Table I: Coupling constants and Gegenbauer coefficients entering the p— meson DAs, at the scale /^o = 1 GeV. Note that 
in Ref. [S^] the normalization are such that /j^^'^ggj — "^p /g^^crc]- 

where 6o = {UN, - 2Nf)/3 , 7„ = 4Cf + 2) + 7^ - | - and ^P{n) = -je + Eltl ■ For 

the f^p coupling constant, the evolution is given by 

with — — ^ +3Cg {Cg — Nc). The couplings and w^Jj i}(m^)/3^(m^) enter a matrix evolution equation 
[53I I. Defining 

v{^^^') = I I (B4) 



it reads 



with given by 



V{n^)=L{iiYi^^'>V{l) (B5) 



^ 3 



7 r< 2 



3^ a 



(B6) 



Hence we get the dependence of f^p and w^^g 1} by diagonalizing the system. In Fig. [9] we display the three 
independent DAs (^1 (left panel), S (center panel), M (right panel). Figure [TU] shows the DAs (^33 (left panel) 
and (pA (right panel) , and Fig. [TT] shows the DAs (left panel) and 93 J (right panel) as a function of their 
longitudinal variables. 

These figures exhibit the non-negligible effect of QCD evolution on DAs, in particular, for the genuine twist-3 
contributions. 



Appendix C: Integration over k for Too and Toi 

In this appendix we describe the method used to evaluate the integrals over k in Eqs. ([551 [5^ when no 
infra-red cutoff is imposed. Let (ui,U2) be the orthonormal basis such as r = r Ui, and then k — r = {ki — r, ^2). 
In that case, a residue method, as described in Ref.ji^, can be applied for the ki integration. In this basis, the 
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91 




Figure 9: The three independent DAs. Left panel: (pi{y; ^J?) as a function of y\ the red (dotted) hne is for the asymptotical 
DA, the blue (solid) line is for /^^ = 1 GeV^, and the blue (dashed) line is for = 25 GeV^. Center panel: S{yi,y2 = 
0.6; fi^) as a function of j/i; the red (solid) line is for fi^ = 1 GeV^ and the blue (dashed) line is for /i^ — 25 GeV^. Right 
panel: M{yi,y2 ~ 0.6; /x^) as a function of yi; the red (solid) line is for — 1 GeV^ and the blue (dashed) line is for 
I? = 25 GeV^ 



2.0 




Figure 10: The two DAs ip^ (left panel) and ipA (right panel) as a function of y. The red (dotted) line is the WW 
contribution with = 1 GeV^, the blue (dash-dotted) line is for the WW contribution with /i^ = 25 GeV^, the red 
(solid) line is the genuine contribution with /^^ = 1 GeV'^, and the blue (dashed) line is the genuine contribution with 
At^ = 25 GeV^. 



amplitudes read 



isCf2QAB f' _ ^ 2^ P^, P^, 1 1 (fci - r/2)^ + fc| - (r/2)^ 



00 = . / / dyyyipi{y,n ) / dki rffca 



(2^)2(Af2 + (r/2)2) J, ™ - ' J_^-^ J_J- kl + kl (fci - r)2 + kl {k, - r/2f + fc| + NP 

11 1 1 



{yr)'^+yyQ^ {yr f+yyQ^ {ki - yrY + k^ + yy Q"^ {kl - yrY + k^ + yy Q"^ 



(CI) 



and 



isCpAB ( 2, r r 1 1 (fci - r/2r + kl - (r/2)2 

^01= ,^ x9/,,f9 , / /n^9^ / dy{y-y)'Pi{v,iJ. ) dki dfca- 



(27r)2(M2 + {r/2y)Jo " "^''"'^ 7-oo J-oo kl + kl (fci - rf + kl (fci - r/2)2 + fcf + 
yrei yr ei ei{ki - yr) + e2k2 ei(fci - yr) + 62/02 



{yrY+yyQ"^ {yr)^ + yyQ^ {ki ~ yr)^ + kl + yy Q"^ {ki - yr)^ + kl + yy ( 



(C2) 



where ei and 62 are the components of the transverse polarization of the 7* in the basis (ui,U2). We define the 
integrands /oo and /oi as 

tsC 2Q AB ^00 /"OO 

^ {2.Y{AP + {v/2Y) X '^yyy^^^y■^^'") j_J^^ j_jk,F,,{kiM.y) m 
= i2.nXlm r - ^) ^^^^^ /I /I ^) ■ ^^^^ 
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We perform a shift of the variables: k — (x + Ui ) | i.e fci — )■ xi , fc2 — >■ a;2 ■ The shift symmetrizes and rescales 
the momenta of the gluons, which are zeros for xi — ±1 and X2 = 0. The integrands then read 



1 



1 



1 xi + x^-i 



r2 {xi + 1)2 + xl (.Ti - 1)2 + xl g{A'P) + xj 



{yrY + J/J/Q2 ^yj-y _^ yyQ2 J.2 



1 



(a;i + (y - y))2 + g{yy Q2) [xi - [y - y)Y + g{yy Q'^) 



(C5) 



and 



fm{xi,X2,y) = 



x\ + x^-l 



r2 (xi + 1)2 + xl (xi - 1)2 + xl g{M^) + x\ 



yrei ^ yrei 



{yrY + yyQ'^ [yrf + yyQ"^ r 



xi + {y - y) 



Xl 



(y-y) 



(xi + {y-y)f + g{yy Q^) (xi ~ [y ~ y)Y + g{yy Q'^) 



(C6) 



with 



4v + r'^xl 



(C7) 



and /oo , /oi being defined such that 



Tf 



00 



isCp^QAB 
(2^)2(M2 + (r/2)2) Jo 



OO /'OO 



dyyyfi{y,Q ) / dxi / dx2 fooixi,x2,y) 



OO ^ — OO 



(C8) 



and 



Tf 



isCpAB 



01 



(27r)2(Af2 + (,/2)2)y^ 



OO /"OO 



OO -/ — OO 



dy (y - y)fiiy,Q ) / rfa;i / da;2 /oi(a;i,a;2,?/) . (C9) 



1. Ratio Toi/Too without any infrared cutoff 

Since the integrands in (jC8|) and (|C9| oscillate quickly, we have used a method where the integration over xi 
can be analytically performed in order to avoid numerical integration issues. We now integrate over the variable 
Xi using the residue method. The poles of the integrands are the same. The poles enclosed in the below contour 
line for X2 > are 



Xl = ±{y - y) - iVgiQ'^yy) i xi^±l- 1x2 , xi = -i^Jg{AP) 



(CIO) 



18 



As the integrands are symmetric under X2 -H- —X2, the resuh is the same for X2 < 0. The remaining integrals 
over X2 and y are then 

Too(r,Q,M) = — / dyyy^i(y,Q^) / dx2h2^7T) } -Resdfoo] (Cll) 



(27r)2(M2 + (r/2)2) 



pi nOO 5 

/ dyyyipi{y,Q^) / dx2i~2in)^ Res^[foo] 

Jo J-oo j^]^ 



To,{r,Q,M)^ (27r)2(M2^+(r/2)2) y„ fa ' Q') rfa;2(-2z^) ^ i?es,[/oi] , (C12) 

where the exphcit expressions for the residues i?esi [/oo] and Resi[f qi] are too lengthy to be displayed here. We 
then integrate numerically over X2 G [0,oo] and y and finally multiply the result due to the symmetry of the 
pole structure by 2. 

2. Ratio Toi/Too with an infrared cutoff 

In order to integrate over ki and k2 with an infrared cutoff A, we will change the variables xi and X2 by the 
distances of the point {xi,X2) from the singularities in (-1,0) and (1,0). Let bi and 62 be these distances such 
that 

bl = {xi + lf+xl hl^ {xi~lf +xl. (C13) 
We have two solutions for [xi, X2), one restricted to the upper half-plane, 

xi^^{hl-bl), X2 = \^{hi +b2 + 2)(6i -f 62 - 2)(6i - 62 + 2)(-6i + 62 + 2) , (C14) 

and one restricted to the lower half-plane, 

xi^\{bl-bl), X2 = ~\yj{hi +b2 + 2){bi +b2- 2)(&i - 62 + 2)(-&i + &2 + 2) . (C15) 



We can restrain the computations to the upper half-plane because the integrands of (|C9p and (jCSp are invariant 
in X2 ^ —X2- The existence of both solutions requires that 

(61 + 62 + 2)(6i + 62 - 2)(&i - 62 + 2)(-6i + 62 + 2) > (C16) 

or, equivalently, at fixed &i, f'2 ^ ^1 + 2 and &2 ^ l^i ~ 2| . This is the condition for the two circles centering in 
(-1,0) and (1,0) and of radius 61 and 62, respectively, to cross each other. The Jacobian of the transformation is 

>0 (C17) 



V(6i + 62 + 2)(6i +b2- 2)(6i -b2 + 2)(-6i +62 + 2 
The infrared cutoff is included by a representation of the Heaviside distribution: 

e(b,-p.^uk)^^^-±^^-^, (CIS) 

where /3cut is the cutoff for bi and 62- The link with the IR cutoff in GeV is A = ficut^- This ensures the stability 
of the numerical evaluation of the 61 and 62 integrations. Then the amplitudes read 

is C 2C^AI3 /"oo /'61+2 

Too(A;fc) = x2/.f2 I t /o^2^ / yy ipi{y; fi'^) dbi rf&2 6'(&i - ^cnt(A); fc) 6'(62 - /3c«i(A); fc) 

(27r)2(M^ -I- (r/2)^) J\b,-2\ 

26162 



v/(6i + 62 + 2)(6i + 62 - 2)(6i - 62 + 2)(-6i + 62 + 2) 



/oo(2/, 2:1(61,62), 2:2(61, 62)) (C19) 
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^oi(^; ^) = (0 w?J2^^fw9^2^ /' (y - v)^^^v^ '^^i /''^' ^^^20(^1 - fc) (C20) 

(27r)2(M2 + (r/2)2) 7„ 7„ 7|j^__2| 

X e{b2- l3cutW;k) ^ ^^^^ ^^ /oi(y,a;i(6i,62),a:2(bi,62)) • 

+ ^2 + 2)(5i + &2 - 2)(6i - 62 + 2)(-&i + 62 + 2) 

The integrations are performed numerically over 61, 62, and y. The constant fc is chosen to be equal to 10. The 
width of model (|C18|) of the Heaviside distribution equals f = jq GeV. This ensures the stability of the 
computation without significantly affecting our results. 
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